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We calculate the short-range exchange-correlation energy of the uniform electron gas with two 
modified electron-electron interactions. While the short-range exchange functionals are calculated 
analytically, Coupled-Cluster and Fermi-hypernetted chain calculations are carried out for the cor- 
relation energy and the results are fitted to an analytical parametrization. These data enable to 
construct the local density approximation corresponding to these modified interactions. 
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I. INTRODUCTION 

In electronic structure calculations using density 
functional theory (DFT) Q in the Kohn-Sham (KS) 
scheme Q , the central quantity that needs to be approxi- 
mated is the exchange-correlation energy functional. The 
vast majority of approximations for this functional orig- 
inates from the local density approximation (LDA) [2] 
consisting in locally transferring the exchange-correlation 
energy of the uniform electron gas to the inhomogeneous 
system of interest. Actually, it is been realized for long 
that the LDA can describe accurately short-range cor- 
relation effects but is inadequate for long-range corre- 
lation effects (see, e.g., 3). This observation lead to 
the development of the first gradient corrected function- 
als d H, B, 0, M, with the basic idea that the long- 
range contribution to the exchange-correlation energy of 
the uniform electron gas must not be transferred to in- 
homogeneous systems. 

Guided by the same idea, it has been proposed [T^, EH, 
El, EH EH Ell to describe only the short-range electronic 
correlations effects by a density functional, and leaving 
the remaining long-range correlations effects to a more 
appropriate method like Configuration Interaction. Con- 
cretely, the method is based on a decomposition of the 
true Coulomb electron-electron interaction as 



long-range interaction 



<(r)+^(r), 



(1) 



where v£ e (r) is a long-range interaction and v£ e (r) is the 
complement short-range interaction. This separation is 
controlled by the parameter /i. In previous studies [HI, 
Ea, E3l, the error function has been used to define the 



<4,crfM 



erf [fir) 



(2) 



referred to as the erf interaction. More recently [16[ , we 
have used a sharper long-range/short-range separation 
with the erfgau interaction 



ee,crfgau 



(r) 



erf(/ir) 2/i 



(3) 



The method then consists in finding the ground-state 
multi-determinantal wave function of a fictitious sys- 
tem containing only the long-range part of the electron- 
electron interaction V£ = J2i<j v ee( r ij) anQl having the 
same density n than the physical system. The total 
ground-state electronic energy of a physical system is 
then given by 



E= (WT + V^ + VneW 



EH c [ 



(4) 



where T is the kinetic energy operator, V ne is the nuclei- 
electron interaction, is the short-range Hartree energy 
and E£ c is the short-range exchange-correlation func- 
tional defined as the difference between the standard 
KS exchange-correlation energy E xc and the long-range 
exchange-correlation energy E£ c associated to the inter- 
action v%„ 



(5) 
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Eq. (HJ) provides an exact decomposition of the total en- 
ergy into a long-range component written in a wave func- 
tion formalism and a remaining short-range component 
expressed as a density functional. In particular, there 
is no double counting of correlation effects. The only 
unknown quantity in this approach is the short-range 
exchange-correlation functional E£ c which is not the 
usual exchange-correlation functional of the KS scheme 
R 



2 



For a reasonable long-range/short-range separation (/j, 
not too small), E£ c essentially describes short-range in- 
teractions, and it is therefore expected to be well ap- 
proximated by the LDA corresponding to the modified 
interaction 



n(r)e^ c {n{r))dr. 



(6) 



In Eq. ([5]), £f c is the short-range exchange-correlation 
energy per particle obtained by difference from the 
exchange-correlation energies per particle of the uniform 
electron gas with the standard Coulomb e xc and with the 
erf or erfgau interaction e l f c 



(7) 



As for the original LDA in the KS scheme with the 
Coulomb interaction, knowledge of ef, c is crucial to apply 
the LDA to the short-range exchange-correlation func- 
tional. In this paper, we give the expressions of this 
short-range exchange-correlation energy per particle of 
the uniform electron gas with the erf and erfgau modified 
interactions. While the exchange part can be calculated 
analytically, the correlation are derived from Coupled- 
Cluster and from Fermi- hypernetted chain calculations. 

Atomic units will be used throughout this work. 
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FIG. 1: Ratio of the erf (solid line) or erfgau (dashed line) 
short-range exchange energy per particle to the exchange en- 
ergy per particle with Coulomb interaction in the uniform 
electron gas ^(r^^^j) with respect to fir s . In order to 
compare the two interactions, a scale factor has been ap- 
plied on the interaction parameter of the erfgau interaction: 
fi — > 3.375/i. 



where E x is the usual KS exchange energy, a„. cr f = 
l/(2n + l) for the erf interaction and a„ j0r f gau = l/(2n + 
1) — 1/3™ for n > 2) for the erfgau interaction. Ex- 
cept for the term linear in /x, the expansion for fir s — ► 
of the erf short-range energy per particle of the uniform 
electron gas does not exhibit the same behavior 



II. SHORT-RANGE EXCHANGE ENERGY 

The short-range exchange energies per particle £ x {r s ) 
of the uniform electron gas for the Wigner-Seitz radius 
r s = (3/(47rn)) 1 / 3 with the erf and erfgau interactions 
are calculated analytically (see Eq. IA9I and IA12I of Ap- 
pendix [A"]) . The inverse of the interaction parameter, 
l//i, represents the range of the modified interaction and 
has to be compared with r s which is the characteristic 
length for exchange. Thus, the relevant variable for the 
exchange energy is actually [ir s . Fig. Q] shows the ra- 
tio of the short-range exchange energy per particle with 
the erf and erfgau interactions to the exchange energy 
per particle with the Coulomb interaction e£(r s )/e x (r s ). 
In order to compare the two interactions, a scale fac- 
tor is applied on the parameter of the erfgau interaction 
fx — > (1 + dty/i) 1 / 2 /i pa 3.375/i so that the erf and erfgau 
exchange energies have the same asymptotic behavior for 
fir s — > oo (see below). 

It has been shown that for a finite system where the 
exchange contribution to the second-order density ma- 
trix n2,a;(ri,r2) decays exponentially with r\2, the short- 
range exchange energy can formally be expanded 
around /i = into an odd series in /i 16] 



1 °° 



7t * — ' n! 

n=0 



( l)"fln 2w+l 



J J n 2 ,z(ri,r 2 )ri2dridr 2 , 



(8) 



r s 1 ( 3 V /3 2 2 . , 
e x {r s ) + — fi - ( ^-j ) r st i + —r sl i 



V 2tt 4 J 



exponential terms 



(9) 



Similarly, the expansion corresponding to the erfgau in- 



teraction is 



£ x, erfgau ( r s) ~ £ x{r s ) + 



2V3 - 3 



r s fi 



2(9 



(187r 4 )!/3 
exponential terms . 



81tt 



(10) 



These different behaviors of the short-range exchange en- 
ergy in the uniform electron gas and in a finite system 
is consistent with the important LDA error arising at 
fi — 0, i.e. for the standard DFT within the Kohn-Sham 
scheme. 

The short-range exchange energy of a finite system can 
also be formally expanded for /i — > oo into the asymptotic 
series |17| 



A 

V ^ (2n)!(2n + 2)/i 2 »+ 2 



n=0 



<2n) 



(r,r)dr, 



(11) 

where n^ 2 ™ (r, r) are the exchange contribution to the 
on-top second-order density matrix and its derivatives, 



*-n,erf 



r( :2 Y^) for the erf interaction and A n el f g 



2 

r(2±2)-3 I1 * 1 r(2±2) + 2x3 il ^r(^) for the erfgau in- 
teraction. The asymptotic expansions of the short-range 
energies per particle of the uniform electron gas for large 



3 



fx do have the same form 



<erf( r «) 



(r ) 

a:, erfgau V sj 



3 1 



16 rf /i 2 



3^ 2 \ 1/3 27 1 



640 rf^i 4 



. (12) 



3(l + 6\/3) 1 
16 rffi 2 

3tt 2 \ 1/3 27(1 + 36V3) 1 



640 rf/x 4 



•(13) 



Again, this is consistent with t he q uality of the local den- 
sity approximation for large \i O E3] ■ 



III. SHORT-RANGE CORRELATION ENERGY 

The long-range correlation energy per particle s^(r s ) 
with the erfgau interaction has been computed for sev- 
eral values of r s (from r s = 0.2 to 10) and fx (from fi = 
to 25). Coupled-Cluster calculations with double exci- 
tations (CCD), according to a method introduced by 
Freeman [13], have been performed (see Appendix iBjl. 
as well as Fermi-hypernetted-chain (FHNC) calculations 
(see Appendix |C|) . Data for the erf interaction is already 
available [11] . 

Once the long-range correlation energy per particle 
e^(r s ) is obtained, the short-range correlation energy per 
particle s^(r s ) is expressed as 



e£(r s ) - e c (r s ) 1 - 



r°°(r s ) 



(14) 



where s c (r s ) is the correlation energy per particle of 
the uniform electron gas with Coulomb interaction taken 
from the usual parametrization of Vosko, Wilk and Nu- 
sair (VWN) [ljf. According to Eq. (H}, e£(r s ) correctly 
reduces to the VWN value for fx = and vanishes for 

fi — * OO. 

The erf and erfgau short-range correlation energies per 
particle with respect to fi for r s = 0.5 and r s = 2 are 
plotted in Fig. [5] and [3] For erfgau, the differences be- 
tween the results from the CCD and FHNC calculations 
are visible only for a high density (r s = 0.5). Surpris- 
ingly, both methods diverge with the erfgau interaction 
when fi-Jrl < 1, explaining the absence of points be- 
tween (1 = and fi w 0.7 for r s — 2, or between fx = 
and // « 1.4 for r s = 0.5 in Fig. [3J We connect this 
behavior to the attractive character of the erfgau inter- 
action for small fi (see Appendix ID]). In practice, the lack 
of accuracy of the LDA correlation functional for small 
fx because of these missing points does not represent a 
serious problem since the LDA exchange functional pro- 
duces large errors anyway in finite systems in this domain 
of fx, as suggested by its incorrect expansion as fx — > 
(Eq.HDJ). 

It has been shown [1J, [l7| that the leading term in the 
expansion of the short-range correlation energy for large 
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FIG. 2: Short-range correlation energy per particle (dots) of 
the uniform electron gas for the erf interaction with respect 
to the interaction parameter fi for r s — 0.5 and r B = 2 com- 
puted with the CCD method. The analytical parametrization 
(Eq. I18[) is represented by a dashed line. 
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FIG. 3: Short-range correlation energy per particle of the 
uniform electron gas for the erfgau interaction with respect to 
the interaction parameter fi for r s = 0.5 and r s = 2, computed 
with the CCD method (dots) and with the FHNC method 
(triangles). The analytical parametrizations (Eq. I18[) using 
the CCD data and the FHNC data are represented by the 
long-dashed and short-dashed lines, respectively. 



fi is 



Ctt 
2f~i 



^2 / n 2:C (r,r)dr + 



(15) 



where n.2, c (r,r) is the on-top correlation pair density for 
the full Coulomb interaction, C = 1 for the erf interac- 
tion and C = (1 + 6\/3) for the erfgau interaction. For 
the uniform electron gas, n2, c (r,r) can be expressed in 
term of the on-top pair-distribution function go(r s ) so 
that the short-range correlation energy per particle has 
the following exact behavior for fi — ► oo 



i / \ 3C / , 1 



(16) 



An estimation of ga{r s ) which includes the correct limits 
for r s — > oo and r s — > was given by Burke, Perdew and 
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FIG. 4: ^ 2 e^(r s ) computed from the CCD method (dots) and 
from the FHNC method (triangles) for the erfgau interaction 
and r 3 = 2. The horizontal line is the exact limit for /i — > oo. 



Ernzerhof 2C 



ff0 (r s ) = ^(( 7 + r s ) 3 / 2 +/?) 



(17) 



with D = 32/(3tt), A = 
7 = 4.7125. Notice that 



3.2581, [3 = 163.44 and 
with this definition, < 
go(r s ) < 1/2. In Fig. 2J we have plotted /x 2 e^(r s ) com- 
puted with the CCD and FHNC methods with respect 
to fi for r s = 2. This plot actually illustrates a general 
trend: for large values of r s , the correlation energy per 
particle computed from the CCD method does not ex- 
hibit the correct behavior for /j, — ► oo. On the contrary, 
the FHNC method seems to perform better in this limit 
in spite of an important numerical noise. However, for 
small r s (r s < 1), the CCD method becomes exact, since 
it reduces to the Random Phase Approximation (See Ap- 
pendix [Bj, and thus respects the /i — > oo limit. 

The short-range correlation energies per particle for 
the erf and erfgau interactions are represented by the 
analytical parametrization 



f(r s ) 



-X r >) 



1 + ci(r s )jj, + c 2 {r s )n 2 ' 



where c\(r s ) is determined by a fit 



ei(r s 



u\r s 



u 2 rt 



1 + v\r s 



(18) 



(19) 



with ui = 1.0271, u 2 = -0.2302, vi = 0.6197 for erf 
«i = 0.3916, u 2 = 0.0223, Ui = 0.9105 for erfgau 
using the CCD data and u\ = 0.4795, u 2 = 1.0094, 
v\ = 10.1247 for erfgau using the FHNC data, and c 2 (r s ) 
is imposed by the exact limit for /x — > oo 



c 2 (r s ) 



8r2e c (r.) 



3C(g (r s ) - 1/2) 



(20) 



These analytical parametrizations for erf and erfgau in- 
teractions are represented in Fig.[^]and[3]for r s = 0.5 and 



2. The two parametrizations for the erfgau interaction 
differ only at small r s where the imposition of the exact 
fi — > oo limit make both expressions close to the CCD 
data. 
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APPENDIX A: EXCHANGE ENERGY OF THE 
UNIFORM ELECTRON GAS WITH MODIFIED 
INTERACTION 

The uniform electron gas can be considered as a sys- 
tem of N electrons in a box of volume Q with a uniform 
background of positive charge to ensure neutrality, stud- 
ied in the thermodynamic limit (i.e N — > oo and — > oo 
such that the density n = N/Q remains constant). This 
system is described by the electronic Hamiltonian 



H = H + Hi, 



(Al) 



where Hq — T is the kinetic energy operator and H{ nt is 
the electron-electron interaction which can be expressed 
by its Fourier expansion 



(A2) 



i<j k^O 



where v e e(k) is the Fourier transform of the (modified) 
electron-electron interaction v ee (r). The term k = 
corresponding to the Hartree energy has been removed 
since it cancels out with the background energy and the 
electron-background interaction energy, provided that 
the same modified interaction has been applied to all 
these terms. The exchange energy corresponds to the 
first-order correction 



E x = (M„t|$> 



(A3) 



where $ is the ground-state wave function of Hamiltonian 
Hq (a Slater determinant of plane- waves). It has been 
shown [2l[ that the exchange energy reduces to 



E r 



i 



(l-x)dq, 
(A4) 

with k F = {^n) 1 / 3 , x = q/(2k F ), 6(y) = 1 if y > and 
0(y)=Qify<O. 

The Fourier transform of the Coulomb interaction is 



4ir 

^ee,coul(Q) 2 s 



(A5) 



so that Eq. (|A4[) leads after integration over x the well 
known exchange energy per particle e x = E x /N 



^■a:,coul(^*s) 



3 / 18\ 1/3 1 



(A6) 



5 



with r s = l/(akp) and a — (4/(97r)) 1 / 3 . For the erf into Eq. (|A4|) leads to the long-range exchange energy 
interaction, inserting the Fourier transform per particle 



An 

<e,erf(?) = ~ e 



-9 2 /(V 



(A7) 



c,crf( r s) 



18 



1/3 



-A Lft&i (J^j + (2A - 4A 5 )e- 1 '^ - 3A + AA 3 ^j , 



where A — /u/(2fcj?). The short-range exchange energy per particle is then 



£x\n) 



"a;, erf 

\>A 1/3 l 



(n) 



^ ^erf + (2A - ^^e" 1 /^ 2 ) - + 



(A8) 



(A9) 



Similarly, the Fourier transform of the erfgau interaction so that the long-range exchange energy per particle is 
writes 



u 

"ee,erigau VI/ ^2 



4?r - 9 2 /(V) _ 6 ^ 7r 



e -3 9 2 /(V) ) (aio) 



'a.erfgauV 4 V \ TT 2 I r 



18\ 1/3 1 



A C^erf f -L) + (24 - 4A 3 ) e - 1 /( 4A2 > - 3A + 4A ; 



-A L/Hed (J^j + (2B - lGB^e- 1 /^ - 6B + 16B 3 ^j 
and the short-range exchange energy per particle is 

- - A ( V^erf (^-]+ (2A - AA 3 ^ 1 /^ -3A + AA 3 



-"X, erfgau 



2.4 



-A ( V^Ferf f±) + (2B - 165 3 )e~ 1/(4s2) - 6B + 16B : 



r 



(All) 



(A12) 



where B = fi/(2^/3kp). 



APPENDIX B: COUPLED-CLUSTER 
CALCULATIONS OF THE UNIFORM 
ELECTRON GAS WITH MODIFIED 
INTERACTION 

For the Coulomb interaction, Freeman has calcu- 
lated the correlation energy of the uniform electron gas 
by summing the ring and screened exchange diagrams 
using the Coupled-Cluster method with double excita- 
tions (CCD). In this appendix, we rapidly give the cor- 



responding equations for an arbitrary electron-electron 
interaction v ee (r). 

The CCD wave function is constructed from the non- 
interacting determinant of plane waves $ through 



*) = e T2 |$), 



(Bl) 



where the excitation operator T 2 is expressed in second 
quantization notation as 



T 2 = *q( k i' k i) a k»+q a k,-q Cl ki Ct k < 



(B2) 



Retaining only the ring diagrams, the amplitudes 



6 



t q (ki,kj) are solutions of the equations (with momen- turn in kp units) 



^q(k^ , kj ) 



V ee (q) 



3ir 2 k F Dq(ki, kj) 



/dk 
_-(t q ( ki! k) + t q (k j( k))0(l - k)9(\k + q| - 1) 



" 18 ^/ / (|^(iq( k - k ) i q( k ,. k ')+i q (k. i ,k')i q (k J ,k)) 



(2tt) 3 y (2tt) e 

x0(i - fc)0(i - fc')#(l k + q| - i)^(l k ' + q| - 1> 



(B3) 



r 



with £) q (ki, kj 



- (q 2 + q. (ki + kj ■ ) ) . Compared to the 



original work of Freeman, the only modification appears 
in the Fourier transform v ee (q) of the arbitrary electron- 
electron interaction v ee (r). Once the amplitudes have 
been computed, the correlation energy per particle can 
be calculated by 



(B4) 



where e c .dir and e c . ex are the direct and exchange contri- 
butions given by 



£c,dir — 187T k F 



f/q 



dk,; 



-W ee (Q)iq(k 4 , kj 



(2tt) 3 J (2tt) 3 J (2tt) 3 
x0(l - h)Q{l - fc,)0(|k 4 + q| - l)0(|k,. + q| - 1), 



(B5) 



dq 



" !, " lA '' ' (2tt) 3 J (2^) 3 
x0(l - ki)6{\ - k 3 )9(\k, 



dki 



dk 

j^Veedki + kj +q|)t q (k i ,k J ) 
-l)0(|k i+ q|-l). 



(B6) 



The direct contribution, corresponding to the ring dia- 
grams, is the usual correlation energy within the Random 
Phase Approximation (RPA) . The exchange contribution 
includes additional screened exchange diagrams. In the 
high-density limit (r s — > 0), the method reduces to the 
RPA and thus becomes exact. 

In practice, it is convenient to introduce the interme- 
diate quantity 

/dk 
^kMl-feMlk + ql-l), (B7) 

and to perform the integration by Gauss-Legendre 
quadrature. Equation (|B3|) is then equivalent to 



B, 



(B8) 



E, Dim 127T 



37T 2 kF kp 



(BIO) 



where — f ee (g)/Z? g (k i ,k J ) and W m are the quadra- 
ture weights. As Bi actually depends on the Tj's coeffi- 
cients, Eq. (|B8[) have to be solved iteratively. 



APPENDIX C: FERMI-HYPERNETTED CHAIN 
THEORY FOR HOMOGENEOUS SYSTEMS 



Similar in spirit to the CCD approach (|Bip , the Fermi- 
hypernetted chain (FHNC) method [12] is based on an 
approximate product ansatz for the wave function 



with Tj — T q (kj) and 



A.; 



k F ^ 



W m D, 



-W jDij , 



(B9) 



*(ri,r 2 ,...,rjv) = exp 



i<j 



$(ri,r 2 ,...,rjv), 
(CI) 



7 



where the correlation factor, called Jastrow factor, acts 
on a single Slater determinant For homogeneous sys- 
tems the pair-correlation function u 2 depends only on 
the inter-electron coordinate r^. The close relation- 
ship between CCD and FHNC methods is not restricted 
to a purely formal analogy between the pair-correlation 
function u 2 and the CCD excitation operator T 2 . This 
topic has been extensively discussed in a review article 
by Bishop 23]. For bosonic systems, both methods are 
actually equivalent on a certain level of approximation. 

It is an important feature of the Jastrow ansatz that 
the exact short- and long-range asymptotic behavior of 
a homogeneous system can be expressed as simple func- 
tions of the inter-electron coordinate. In the case of a 
Coulomb potential, Kato's cusp condition for electrons 
with antiparallel spin imposes a constraint on the first 
derivative of the pair-correlation function 



du 2 {ri 2 ) 



dr 12 



ri 2 =0 



1 

2' 



(C2) 



which can be exactly represented by a Jastrow factor. 
We discuss below how the modified interaction affects 
the short-range behavior of the Jastrow factor. The long- 
range asymptotic behavior of electron correlations is well 
described by the RPA approximation [24[. It provides 
an explicit asymptotic expression for the pair-correlation 
function 



lim u 2 {r 12 ) = : 

ri2-»oo w pl r 12 



(C3) 



where the plasmon frequency w p j = \J Aim of the elec- 
tron gas enters into the denominator. This asymptotic 
behavior can be reproduced by the FHNC method [25| . 

For a given pair-correlation function, the FHNC equa- 
tions represent a nonlinear system of equations between 
"nodal" N(r 12 ), "non-nodal" X(r i2 ) and "elementary" 
E(r 12) functions. Each of these functions can be ex- 
pressed as an infinite sum of certain types of diagrams 
build up from the pair-correlation function and the one- 
particle density matrix of the noninteracting system. 
Some of these equations are conveniently expressed in 
coordinate space, the others in momentum space. The 
system of equations is underdetermined and requires an 
a priori knowledge of the "elementary" diagrams in order 
to get a unique solution. In a series of papers Krotscheck 
developed a consistent approximation scheme for the 
FHNC equations [26|, [27l |28||. which preserves the cor- 
rect asymptotic behavior on each level of approximation. 
We have used the FHNC//0 method which corresponds 
to the lowest level of approximation, where "elementary" 
diagrams are neglected altogether. The FHNC//0 equa- 
tions are given by 



T d d{ri 2 ) ■= X dd (r 12 ) + N dd (r 12 ) 

= exp[2u 2 (r 12 )+iV dd (r 12 )]-l, (C4) 



N dd (k)=X dd {k)S F {k)T dd {k), 



where S F is the liquid structure function of the noninter- 
acting system. We have used the dimensionless Fourier 
transform 



f(k) — n J drf(r) exp(ik • r) 



(C6) 



The "nodal" and "non-nodal" functions provide a link 
between the Jastrow factor and the liquid structure func- 
tion of the interacting system 



S(k) = S F (k) + S F (k) 2 f dd (k), 



(C7) 



which is essentially the Fourier transform of the pair- 
density n 2 (ri 2 ) 

S(k) = 1+n J dr 12 (n 2 (ri 2 )/n 2 - l) exp(ik-r 12 ). (C8) 

This connection enables an approximate variational 
treatment of the Jastrow factor within FHNC theory. 

In the following we want to give a brief outline of 
the FHNC//0 optimization cycles following essentially 
Krotscheck's paper [28| . Starting point is an effective 
particle-hole potential 

V ph (r 12 ) = [1 + r dd (n 2 )] « ee (ri 2 ) + |v[i + r dd (r 12 )] 1/2 

+r d£i (r 12 )a; 7 (r 12 ), (C9) 

which depends, beside diagrammatic contributions, on 
the bare (modified) interaction potential v ee and an in- 
duced interaction u>i. In momentum space the induced 
interaction 



k 2 



1 



S(k) 



S F (k) 



1 



1 



-1 2 



S(k) S F (k) 



, (CIO) 



can be expressed in terms of the liquid structure functions 
of the interacting and noninteracting system. Within 
the high density regime, v ee can be taken as an initial 
guess for V p h- Performing FHNC calculations at suc- 
cessively lower densities it is possible to reach the low 
density regime by taking V p h from a slightly higher den- 
sity as an initial guess in the optimization process. The 
particle-hole potential is related to the liquid structure 
function 

S(k) = - _ (oil) 



l + (A/k^)S 2 F (k)V ph (k) 



In the first step of the optimization cycle Eq. (jCllj) 
is used to get an improved approximation of the liq- 
uid structure function. Using Eqs. (|C7p and (|C10|) it 
is now possible to obtain improved approximations for 
the induced interaction luj and the diagrammatic quan- 
tity Y dd (ri2). These can be used in the second step to 
calculate an improved approximation of the particle-hole 
potential V p h (Eq. IC9|) . The two steps provide a self- 
consistent optimization cycle, which can be repeated un- 
(C5) til convergence has been achieved. Finally we have used 
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FIG. 6: Fourier transforms of the Coulomb interaction (dot- 
ted line), erf interaction (dashed line) and erfgau interaction 
(solid line), each plotted with a interaction parameter /i = 0.5. 



FIG. 5: Jastrow factors at r s = 1 for Coulomb interaction 
(solid line), and erfgau interactions at interaction parameters 
\x = 20 (dashed line), /j, = 10 (dotted line) and fj, = 5 (dotted- 
dashed line). 



the FHNC Eqs. (|C4j and (JC5J) in order to obtain the 
optimized FHNC//0 Jastrow factor. 

Jastrow factors for the Coulomb interaction and the 
long-range erfgau interaction are shown in Fig. [5l With 
decreasing value of the interaction strength /i, the short- 
range part of the Jastrow factors is modified; it changes 
from a cusp at r%2 = for \x — > oo into a smooth behavior 
for any finite fi. For small values of fx, a local minimum 
appears at an intermediate distance. As expected, the 
long-range behavior of the Jastrow factor is not affected 
by the modified interaction. 



APPENDIX D: DIVERGENCE OF 
CALCULATIONS ON THE UNIFORM 
ELECTRON GAS WITH MODIFIED 
INTERACTION 

With the erfgau interaction, CCD and FHNC calcula- 
tions of the uniform electron gas diverge for small values 
of fi and r s . This is due to the particular form of the er- 
fgau interaction. In fact, whereas the Fourier transform 
of the Coulomb or erf interaction is always positive, the 
Fourier transform of the erfgau interaction 



„ . (n) = !V 9 2 /(V) - 

u ee,crtg I au\ L l I o ° o 

q z [l 1 



-3<2 2 /(V 



(Dl) 



can be negative (see Fig. [6|. For small n, the negative 
part of f e e,erfgau(<z) is not negligible, introducing an at- 
tractive contribution to the electron-electron interaction. 

It is possible to estimate the domain of /i and r s for 
which CCD calculations do not converge because of this 
attractive interaction. Let's consider the "state-average" 
approximation to the CCD equations proposed by Bishop 
and Liihrmann [29j . This model consists in neglecting 



the exchange contribution to the correlation energy and 
averaging over the occupied momentum and the 
equations given in Appendix [B] The correlation energy 
per particle of the uniform electron gas is then written 
as (with momentum in kp units) 



4^ 



dqq 2 v ee (q)P(q)(t q 



(D2) 



(i 



where P(q) = 3q/4 - q 3 /16 if q < 2, P(q) = 1 if 
q > 2, v ee (q) is the Fourier transform of the electron- 
electron interaction and (t q ) is the average of the ampli- 
tude t q (ki,k,) which is solution of the equation 



(t g ) = 



1 



iix^kp 



v 1 ) 



(D3) 



where (D q 1 



is the average of the inverse of D q (ki,kj 



introduced in Appendix [Bl 

The general solution of (|D3|) is 



(t q ) = 



1 - A(q) + 2A(q) 
A(q)P(q) 



(D4) 



with A(q) = (2(D- 1 )v ee (q)P(q))/(3TT 2 kp). Using the ad- 
ditional approximation (-D" 1 ) ~ (Dq)^ 1 = —P{q)/q 2 , 
one sees immediately that this solution breaks down 
(more precisely, becomes imaginary) if 



MX 



3Tr 2 k F q 2 
' AP{q) 2 



(D5) 



Let's evaluate this inequality in the worst situation where 
v ee (q) and P(q) take their minimum values. The inter- 
action v ee {q) reaches its minimum v m i n m —5.6k F /fi 2 
for q fJ-/kp, and in this domain P{q) < 3g/4 so that 
condition (|D5[) roughly gives 



(D6) 



where r s = l/(akp) with a = (4/(97r)) 1 / 3 has been used. 
In Fig. [71 we have reported the values of \x and r s at 
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U 3 
2 



the limit of convergence for the calculation of the corre- 
lation energy of the uniform electron gas with the erfgau 
interaction, together with the divergence condition (|D6jl . 
Obviously, the domain of divergence is well approximated 
by this condition. 



o 1 — ■ — ■ — ■ — ' — ■ — ■ — ■ — ' — ■ — ■ — ■ — ' — ■ — ■ — ■ — ' — ■ — ■ — ■ — J 

2 4 6 8 10 

r s 

FIG. 7: Points at the limit of convergence for calculation of 
the correlation energy of the uniform electron gas with erfgau 
interaction. The domain of divergence is well approximated 
by [lyfrl < 1 (solid line). 



The divergence of both CCD and FHNC methods for 
these values of \x and r s where the attractive part of 
electron-electron interaction becomes important is rem- 
iniscent of the situation happening in a superconductor 
where the ordinary perturbation expansion breaks down 
for the superconducting phase. 
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